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The motion of granular material in a ball mill is investigated using molecular dynam- 
ics simulations in two dimensions. In agreement with experimental observations by 
RothkcgelEJ we find that local stresses - and hence the comminution efficiency - are 
maximal close to the bottom of the container. This effect will be explained using anal- 
ysis of statistics of force chains in the material. 



1 Introduction 

Milling is one of the most important techniques in mechanical engineering and much 
effort has been done to optimize the comminution efficiency of milling techniques 
with respect to energy, space and time consumptionH. 

Mainly because of their simple construction and application ball milling is a 
wide spread milling technology, particularly in mining. Ball mills are (usually large) 
cylindrical devices of length and diameter up to several meters which revolve along 
their axis. The material to be comminuted moves inside the cylinder. There are 
two different methods: autogeneous comminution where the material is pure and 
heterogeneous comminution where the material is mixed with heavy spheres from 
steel of typical diameter of several centimeters to increase efficiency. Throughout 
the paper we will deal with autogeneous comminution. Unfortunately the efficiency 
of ball mills is not very high and, therefore, engineers did much scientific work to 
increase their efficiency. There exists much experimental knowledge on the operation 
mechanisms of ball mills and on the comminution of granular matter in these mills. 
An overview can be found in ref.a. 

There are still effects, however, which have not been understood yet. The 
present paper deals with the spatial distribution of stress in the material. At first 
glance one could assume that the largest stresses will be observed close to the 
surface of the material where the granular particles have the largest values of relative 
collision velocity. Hence one could assume that large part of particle comminution 
occurs close to the material surface. To increase efhciency of the mill one would 
have to tune the rotation velocity so that the average collision velocity becomes 
maximum. 

Experimental investigations, however, show opposite results: Jlothkegel and 
Rolf measured directly the spatial distribution of intensive impactaj in a ball mill. 
In an experiment using a small almost two-dimensional mill they used instrumented 
balls which flash whenever the acting force at a particle contact is larger than a 
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threshold. The statistics of the spatial distribution of the flashes led to a surprising 
result: the spatial distribution of intensive particle contacts which in realistic ball 
mills might lead to particle comminution has its maximum deep inside the material 
close to the bottom of the cylinder. This result was new and there was no satisfying 
explanation for the measured stress distribution yet. 

In the following section we present the results of a two dimensional molecular 
dynamics simulation of a ball mill. It will be shown that the spatial stress dis- 
tribution is closely related to the properties of force chains in the material. The 
occurrence of such force chains ,^nd even of entire networks of such chains has been 
reported before, see e.g. refs.0H3. 

2 Molecular dynamics model 

The rotation velocity of ball mills is high enough to keep the granular particles in 
the continuous flow regime (for explanation of the regimes seed) and on top of the 
material grains are thrown through the air caused by intensive forcing due to rapid 
revolution of the mill. Therefore static friction does not play an important role for 
the dynamics of the material and, hence, we can apply a sphere model in the MD 
simulations. We just want to remark that in other cases when static behavior of 
granular material becomes essential for the dynamics one has to apply more compli- 
cated grain models of non-spherical shape. For detailed discussions of tiie-differences 
in the dynamics of spherical and non-spherical objects see referenceaHJllj. 

For our simulation we apply a sphericaLMD model for granular materials by 
Cundall and StracktSI and Haff and Wernertil which was used by many authors in 
simulations of rapidly moving granular material. 

Two particles i and j with radii Ri and Rj, positions r^, fj and velocities fi 
and fj feel a force 

F^.^F.^rt + F^'^t (1) 

only when contacting, i.e. if the condition \fi — fjl < Ri + Rj holds. The force 
consist of a component Fjj acting in normal direction n 

F,'J = Y-{R, + R, - \n - I) - mf ■ 7^ • (r^ - f,) ■ n (2) 

and a component F^^ acting in tangential direction ^=(1*^ 0)'^ 

(4) 
(5) 

(6) 

rrii and uji are the mass and the rotation velocity of the i-th particle. Y — 
8T0® g is the Young Modulus, jn = 800 and 7t ~ 3000 are the damping 
coefRcients in normal and tangential direction, and = 0.5 is the Coulomb friction 



F' 



sign [ylfj ■ min ym^f jt 
with 

vlf = (7=*, -fj)- t + R^-LO^ + Rj 

eff rrii ■ rrij 

-' rrii + m.i 
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coefRcient. The parameters have been chosen to give reahstic resuhs in comparison 
of the simulation resuhs with the behavior of a typical granular material. 

Eq. (|^) describes the relative velocity of the particle surfaces at the contact 
point. This velocity consists of the tangential part of the relative velocity of the 
particles and of a term which originates from rotation of the grains. stands for 
the effective mass. 




Figure 1: Snapshots of the simulation for different velocity. The color codes for the local pressure 
Pi (cf. eq. (M)). Black lines connecting adjacent particles mark automatically detected force chains. 
The ball mill revolves with angular velocity Qj = 2Hz (top), Qjj = lOHz (middle), Qjn = 19Hz 

(bottom). 

The Coulomb friction law is taken into account by equation (||) saying that two 
grains slide on each other if the tangential force is larger than fi times the normal 
force. Otherwise they roll mainly. For the detailed discussion oL-this and other 
models for molecular dynamics simulation of granular material seeli3. 

Based on the force (P Newton's equations of motion have been integrated using 
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a Gear predictor-corrector scheme of fifth order, e.g.Q. This method was used in 
many simulations of granular material and has been proven to be very numerically 
stable. 

3 Results of the simulation 

The model described in the previous section was used to simulate a two-dimensional 
ball mill of diameter D = 8 cm filled with — 800 spherical grains. The radii of 
the grains are distributed in the interval [0.05, 0.11] cm. To study the influence of 
the rotation velocity of the cylinder we performed simulations for three different 
velocities ftj = 2Hz = 2/(2 • tt) revolutions per second, fljj = lOHz = 10/(2 • 
7r)revolutions per second, ilm — 19Hz = 19/(2 ■ 7r)revolutions per second. 

At ^ il] — 2Hz the flow at the surface is already continuous and the surface 
is characterised by an inclined line (see fig. |l|, top). For fl = fljj = lOHz the 
material surface is still almost compact and consists of two sections a steep and a 
flat one (fig. |l|, middle). This velocity regime is experimentally known to have the 
best comminution efRciencya. Large rotation velocity SI = fliii = 19Hz leads to 
free flight of the grains at the surface (flg. |^, top) . 

Fig. shows snapshots of simulations with TV — 800 grains for the three men- 
tioned rotation velocities fi/, flu, SI///. The color codes for the local pressure Pi 
acting on the i-th particle 

The index j runs over all neighbors of the ith particle. Grains feeling a pressure 
Pi > lOOOg cm are drawn bold lined. Obviously most of the hardly stressed 
particles are located near the walls inside the bulk oL material, which agrees well 
with the experimental findings by Rothkegel and RolftJ. 

The hypothesis we want to discuss in the following is, that the self organized 
formation of force chains is responsible for the spatial stress distribution and, hence, 
the relevant physical reason for comminution processes in ball mills. Without the 
existence of such force chains, ball mills would not work or at least would work 
much less efflcient. In the following we will discuss the properties of these force 
chains. 

Force chains are defined by a set of simple conditions: Particles are called 
members of the same force chain if: 

1. they feel a pressure (eq. (j^)) of more than lOOOg cm s^^, 

2. each member particle of the chain touches at least one other member, and 

3. for all members k of the chain having two neighbors i and j, which belong to 
the force chain as well, the condition 



rk - f-i rk - r,- 



\rk-ri\ \rk-rj\ 
holds, i.e. these three particles he almost on a line. 



< 0.85 , (8) 
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These three conditions can be checked using a computer algorithm. In fig. |l| 
neighboring particles have been connected by a black line provided that either of 
them belong to the same force chain. Obviously most of the hard stressed (and 
therefore bold drawn) particles are members of force chains. We conclude that an 
essential part of the static and dynamic pressure in the ball mill propagates along 
force chains. Particles which belong to a force chain feel pressure which are up to 
100 times larger than the local average pressure 



1 ^ 

p(av) _ i p 
k=l 



(9) 



where the index k runs over all K neighbors of the particle i which are within a 
certain neighborhood jri — rfc| < 4 • i? and which do not belong to any force chain. 

The occurrence of force chains in granular systems is not_n£w. They have 
been obseri^fid before in other granular systems in experimentsuTS and computer 
simulationgjfi. It can be shown that one of the (at least) three phases of an idealized 
granular system is characterized by the occurrence of force chainsB. 

Fig. H (left) shows the frequency of force chains of lengths L where L is the 
number of particles which join the chain. The frequency decreases almost exponen- 
tially with increasing length. For the highest rotation velocity flm (fig. H, bottom) 
the frequency distribution breaks down for longer chains, which is mostly due to 
less compact bulk of material and finite size of the system. 
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Figure 2: Left: Frequency of force chains of length L. Right: Average of the maximal pressure 
inside one force chain as a function of the length L of the force chain. Top: Qj = 2Hz, middle: 

Qjj = lOHz, bottom: Qm = 19Hz. 



The right hand side of fig. g shows the correlation of the length of a force 
chain and the pressure which acts on the particles which belong to this chain. Of 
particular interest is not the average but the maximum pressure inside a chain. The 
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maximum pressure is criterion whether a grain wiU be comminuted in a miU. As 
soon as one of the particles of a chain has broken the load discharges and the force 
chain vanishes and reappears elsewhere. Therefore here we are mainly interested in 
the maximum pressure. For low velocities fi/ and fi// (fig. ^, top and middle) the 
average value of the maximal force in a force chain of length L grows monotonously 
with the length, i.e. a force chain acts similar as a beam in a framework: all the 
weight and momenta from particles above the chain are supported by the chain and 
propagated through the chain. Therefore one finds in most cases that the lower the 
position of a grain the higher the load. 

In contrast, the corresponding plot (fig. ||, bottom) for high velocity flm shows 
no characteristic dependence of the maximal force on the length of the force chain, 
i.e. for this case the force chains do not act as described in the previous paragraph. 
Thus we conclude that the mechanism of stress propagation in force chains is much 
less efficient in the case of high rotation velocity. 

In fig. I the average number of particles per force chain which feel a pressure 
Pi larger than a given threshold P is plotted versus the length of the force chain. 
This number is directly-related to the breaking probability which can be described 
using WeibuU statisticalH. Again one can see that longer force chains lead to a higher 
probability of comminution. For the highest rotation velocity ^Iju the distributions 
are significantly lower again. 




Figure 3: Average number of grains per force chain, which feel a local pressure Pi larger than a 
given threshold P. Top: Qj = 2 Hz, middle: O/j = 10 Hz, bottom: Qui = 19 Hz. The longer the 
chain the more particles feel pressure Pi > P and, hence, the higher the comminution probability. 
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The time averaged spatial distributions of the pressure Pi for flj = 2 Hz are 
shown in fig. ^ (left part of figure). The color codes for the local pressure Pi, where 
blue means low pressure and red color stands for high pressure. The upper plots 
include only particles which do not belong to a force chain, whereas the lower figures 
include only particles which do belong to a force chain. 




Figure 4: Left: Time averaged spatial distribution of the local pressure for = 2 Hz. The color 
codes for the pressure Pi, blue color means low pressure and red color stands for high pressure. The 
upper figure includes only particles which do not belong to a force chain, the lower figure includes 
only particles which do belong to a force chain. Right: Time averaged spatial distribution of the 
number of grains feeling a pressure Pi > 3000 cm g s~^. The color code for the right figures is 
given by the color bar. Again in the upper figure only grains which do not belong to a force chain 
have been considered, whereas in the lower figure only grains which belong to a force chain join 

the averaging procedure. 

The right hand side of figure ^ shows the time averaged spatial distribution 
of the number of particles, which feel a pressure Pi > 3000 g cm nSec^^. These 
results compare directly to the experiment by Rolf and RothkegelH (see above). 
The maxima of the spatial distribution are located at the bottom of the cylinder 
near the wall which coincides with the results of the experiment. 

Figures || and |6| are equivalent to fig. ^ where the angular velocity of the mill is 
fi// — 10 Hz and flju — 19 Hz, respectively. The color scaling is the same for all 
left side figures and the color code for the right figures is given by the color bar. 

From these figures we conclude the surprising result that direct impacts of par- 
ticles with high relative velocities at or close to the surface do not lead to high 
stresses. Therefore they are almost neglectable for the comminution process. Re- 
gions of high pressure can be found near the walls deep inside the material. Particles 
which do not belong to a force chain feel everywhere in the material a low pressure, 
only particle which are a member of a force chain feel larger stresses, i.e. only par- 
ticles which belong to a force chain have a high probability to break. Fig. |6| shows 
that the mechanism of force chains is of low efficiency for rather high rotation ve- 
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locities. The pressure is much lower than in the other figures ^and which leads 
to a diminution in the breaking rate as observed experimentallyH. 

For O = fliii (fig. ^ the rates resulting from impacts of particles at the surface 
overcomes the rates resulting from the force chains. Nevertheless the rates are 
clearly lower than for slower spinning mill. 

4 Conclusion 

Using two-dimensional molecular dynamics simulation of spherical granular parti- 
cles we investigated the motion of granular material moving in a ball mill which 
spins with different angular velocity j^The numerical results which agree with exper- 
imental data by Rolf and RothkegelEJ lead to an explanation for the experimentally 
observed spatial distribution of pressure in such a ball mill. A network of force 
chains which permanently changes its structure has been observed in the rotating 
cylinder. Detailed and separate analysis of the statistical properties of these chains 
and of the pressure acting on particles which either do belong to a chain, or which 
do not belong to a chain lead to the conclusion that force chains play an important 
role for comminution processes in granular materials. Force chains lead to concen- 
tration of the dynamic and static load of a large amounts of grains in the bulk of the 
material to few grains close to the bottom of the container. These particles finally 
feel very large pressure. The described mechanism of the action oLforce chains is 
able to interpretate experimental data while other explanations failEI. 
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